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By using virial theorem, Helmholtz and Kelvin showed that the contraction of a bound self- 
gravitating system must be accompanied by release of radiation energy irrespective of the details 
of the contraction process. This happens because the total Newtonian energy of the system En 
(and not just the Newtonian gravitational potential energy E^ ) decreases for such contraction. In 
the era of General Relativity (GR) too, it is justifiably believed that gravitational contraction must 
release radiation energy. However no GR version of (Newtonian) Helmholtz- Kelvin (HK) process 
has ever been derived. Here, for the first time, we derive the GR version of the appropriate virial 
theorem and Helmholtz Kelvin mechanism by simply equating the well known expressions for the 
gravitational mass and the Inertial Mass of a spherically symmetric static fluid. Simultaneously, 
we show that the GR counterparts of global "internal energy" , "gravitational potential energy" and 
"binding energy" are actually different from what have been used so far. Existence of this GR HK 
process asserts that, in Einstein gravity too, gravitational collapse must be accompanied by emission 
of radiation irrespective of the details of the collapse process. 
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INTRODUCTION 

It is generally believed that General Relativistic (GR) 
gravitational collapse must be accompanied by the emis- 
sion of radiation such as photons and neutrinos. There 
have been many studies on radiative spherical GR col- 
lapse and we may recall here only a few of them[l]. 
Though all such studies vastly differ in their details, nev- 
ertheless, all of them indicate that the fluid becomes hot- 
ter during the collapse while net amount of radiated en- 
ergy steadily increases. Surprisingly, though, such effects 
seem obvious at first sight, one would wonder whether 
there is a fundamental reason for occurrence of such phys- 
ical effects during gravitational collapse. Put another 
way, whether, irrespective of solution of actual collapse 
equations with their associated assumptions and simpli- 
fications, one can demand from a general perspective, 
that actual gravitational collapse must be accompanied 
by both emission of radiation and heating up of the fluid. 
From thermodynamics perspective, one can ask whether 
the phenomenon of occurrence "negative specific heat" , 
known for Newtonian gravity for a long time, must be 
valid in Einstein gravity too. To appreciate this, let us re- 
call that the specific heat is defined through C — dQ j dT 
where dQ is the amount of heat injected into the system 
and dT is the corresponding increment of temperature. 
Gravitational compression raises the temperature so that 
dT > 0. A negative C would then demand dQ < and 
vice-versa. Hence a negative dQ means loss of heat (ra- 
diation) from the system and vice-versa. 

We emphasize here that, in a strict sense, this phe- 
nomenon of "negative specific heat" is known only for 
weak Newtonian gravity. Intuitively such an effect is 
expected to be more pronounced for a fluid subject to 



much stronger Einstein gravity. But the actual fact is 
that, there is no proper GR theorem which can assert 
that global Einstein gravity too is marked by the same 
"negative specific heat" . 

To highlight this, in Sec 1, we will first review the 
case in the Newtonian gravitation. This would show why 
Newtonian collapse must be accompanied by radiation 
howsoever small it may be and why global gravitation is 
characterized by "negative specific heat" in Newtonian 
case. Then it would be emphasized that a corresponding 
GR derivation is non existent and accordingly we shall 
present an exact GR counterpart of this Newtonian pro- 
cess. We would then automatically arrive at global defini- 
tions GR Self-Gravitational energy and Binding Energy 
from the perspective of global energy conservation of a 
static spherically symmetric fluid. The entire exercise 
will show why, irrespective of the details, GR collapse 
must be accompanied by emission of radiation and an 
increasing fluid temperature. 



NEWTONIAN GRAVITATIONAL COLLAPSE 

As per Chandrasekhar[2], von Helmholtz first proposed 
in 1854 that contraction of self-gravitating bodies should 
emit radiation. Few years later, in 1861, Kelvin[3] elab- 
orated on this process of energy generation which may 
be called as Helmholtz - Kelvin (HK) process. Without 
going into further historical account, the basic physics 
behind the H-K process is reviewed below from a rela- 
tively modern perspective [4, 5]: 

If we consider a spherically symmetrical static isotropic 
fluid in hydrodynamical equilibrium, it follows that 
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E" + 3 I P dV = 47ri? a p b 



(i) 



where, Eg is the Newtonian gravitational potential en- 
ergy, p = p{r) is the isotropic pressure, pb is the pres- 
sure at the boundary r — R, and the volume element 
dV = 4irr 2 dr. For a laboratory gas sphere, it is possible 
to have p(r) w uniform p b . In such a case, Eq.(l) 
would reduce to 

E% + AttR 3 Pb w 4^i? 3 p b (2) 

and which would express the obvious fact that for a lab- 
oratory gas Eg w 0. Note that it is necessary to have 
Pb > in such a case. 

In contrast, if the fluid is assumed to be self-contained, 
i.e., bound by its own gravity then one expects 

Pb = (3) 

and obtains the better known form of Eq.(l): 



E a 



3 pdV = 



(4) 



This is known as static and scalar virial theorem. If the 
adiabatic index of the fluid is 7, then 



P= (7- l)e 



(5) 



where e is the internal energy density. For some ideal 
fluids, 7 may also be considered as the "ratio of specific 
heats" . For example, for a monoatomic ideal gas having 
an equation of state (EOS) p = nkT, where n is num- 
ber density of the monoatomic molecules and k is the 
Boltzmann constant 7 is the ratio of specific heats with 
a unique value 7 = 5/3. For an ultrarclativistic gas with 
particle momenta — > 00 or for a pure photon gas, 7 is 
again the ratio of specific heats having the unique value 
7 = 4/3 [2]. The internal energy for the entire fluid is 



U 



dV 



(6) 



Using Eqs.(5) and (6) in (4), and assuming 7 to be uni- 
form, we have 



E% + 3(7-1) U = 



so that 



U = 



-1 



-E 



N 



3(7-1) 9 
The total Newtonian energy of the fluid is 

E N = E? + U=^- T jE? 



(7) 
(8) 

(9) 



Note that for the above mentioned ideal gas having an 
EOS p — nkT, Eq.(7) reduces to the more familiar form 
E™ + 2U = because here 7 = 5/3. However Eqs.(l-9) 
are valid for any fluid obeying relation (5) and not nec- 
essarily by an "ideal gas" alone, as long as the fluid may 
be assumed to obey an EOS of the form (5). Thus, in 
principle, 7 is arbitrary here subject to general thermo- 
dynamical constraints. 

In Eq.(9), En is the Newtonian binding energy of the 
fluid and must be negative for a fluid which is already as- 
sumed to be bound. For attractive gravity, for any case, 
bound or unbound, one must have Eg < 0. In the present 
case, as soon as we set pb = 0, we imply the system to 
be self-bound, and one must have En < 0. Thus from 
Eq.(9), it transpires that, one must have 7 > 4/3. A 
limiting case of 7 = 4/3 would signify a transition to 
unbound systems. For "unbound systems" one expects 
to have, En > 0. If the system would be unbound, one 
would have pb > and further the system would not be 
in hydrostatic equilibrium. In such a case, one needs to 
use dynamical form of virial theorem to study it. Also, 
by noting Eq.(7), one might think that for unbound sys- 
tems, one might have 7 < 4/3. But this would be an 
incorrect conclusion, because as explained above, for un- 
bound systems, Eq.(7) would cease to be valid. It may 
be borne in mind that 7 is an inherent thermodynam- 
ical parameter and cannot be dictated by gross global 
hydrodynamical behavior of the fluid. 

A limiting value of 7 = 4/3 corresponds to a situation 
when the momenta of the constituent particles of the 
fluid, P = 00, and thus, cannot be strictly realized ex- 
cept for singular situations [6]. For instance note that the 
critical ultrarelativistic White Dwarf of Chandrasekhar 
has R = because it strictly corresponds to a fluid hav- 
ing 7 = 4/3 [2]. 

If, additionally, the fluid obeys a polytropic equation 
of state 



p = Kpi' 



(10) 



where, K and 7 P are uniform over the fluid, it follows 
that 



E 



N 



i GM 2 



9 5-n R 



(11) 



where 7 p = 1 + 1/n and G is the gravitational constant. 
Note that, in general, 7 ^ j p and only if the fluid is 
considered to undergo adiabatic change, one would have 
7 = 7 P . Further, as we would see, all contraction pro- 
cesses are expected to be accompanied by emission of ra- 
diation, and they must be non-adiabatic in a strict sense. 
It may be also mentioned that, in the Newtonian case, 
the fluid density appearing in Eq.(10) essentially means 
rest mass density: p = po- 

If one differentiates Eq.(8), for slow contraction, one 



3 



will have 



dU 
dt 



-1 dE, 



N 



3(7-1) dt 
Also, from Eq.(ll), we see that 



dE, 



N 



dt 



3 GM 2 
5-n R 2 



R 



(12) 



(13) 



Since 7? < for contraction, while the value of E, 
creases during contraction its absolute value 



» de- 
E»\ in- 

creases. From Eq.(12) , we find that, as \E™\ increases 
during such contraction, so does U. However, the amount 
of total gravitational energy released by the contraction, 
\dEg\, is not fully accounted for by the gain in the value 
of U: 



dU = 



3(^1)1^1 < 



N I 



(14) 



For overall energy conservation, it is therefore necessary 
that the rest of the energy gain 



1 



1 



3(7-1) 



\ dE ^ = it^T)\ dE ^ = dQ ^ 



is radiated away by the system. This could have been 
found directly by differentiating Eq.(9): 



dE N 37-4 dE, 



N 



dt 



3(7 - 1) dt 



Using Eq.(13) into above Eq., we find that 



GM 2 



3(7-1) 5-n R 2 



dE N _ 37-4 
dt 



R<0 



(16) 



(17) 



In case of a gas confined in a laboratory by physical 
inclosure, pb > 0. One can also imagine the physical 
inclosure to be a perfect insulator and one may con- 
ceive of a radiationless adiabatic contraction for arbi- 
trary 7 > 4/3. But in an astrophysical context, there 
is neither any physical inclosure nor any perfect insulat- 
ing surrounding. Hence, it appears from Eq.(15) that a 
strictly adiabatic contraction (dQ = 0) would be possi- 
ble only for the idealized case of 7 = 4/3. And Eq.(9) 
would show that, in such a case, one would already have 
En = 0, i.e., they system would be unbound. In reality, 
one has 7 = 4/3 only for pure radiation or when the en- 
ergy of the particles per unit rest mass E* = 00, which 
is possible only for a singular situation in case the "gas" 
is not already a pure radiation [6]. 

Thus, Eq.(17) shows that the total (Newtonian) energy 
of the system decreases for contraction and it could be 
so only if the system radiates appropriate amount of en- 
ergy. Since U increases, the fluid become hotter while it 
radiates (dQ < 0). Therefore a self-gravitating fluid has 
a negative specific heat and this fact is well known. Note 



that this result follows from the Newtonian HK process 
and does not depend on the details of either the physical 
properties of the fluid or the collapse process. 

Since the above result is of generic nature, it is ex- 
pected to be qualitatively valid even in case of strong 
gravity. In fact, even after the introduction of General 
Relativity (GR) into astrophysics, the idea that gravi- 
tational contraction must result into radiation output is 
naturally and justifiably used[l]. While considering the 
configuration of static fluid spheres in GR, Buchdahl[7] 
posed the question whether the amount total radiation 
emitted by gravitational contraction can even exceed the 
initial value of E = Mc 2 itself. The answer was in the 
negative. Yet, for Einstein gravity no exact counterpart 
of Eqs. (12-17) exists. Thus, we cannot assert, as a prin- 
ciple, that self-gravitating matter has "negative specific 
heat" in Einstein gravity too. And we want to address 
this precise aspect in the present paper. 



STATIC FLUIDS IN GENERAL RELATIVITY 

Let us consider a static self-gravitating fluid sphere 
described by the metric 



ds 2 



A 2 (r)dt 2 - B 2 {r)dr 2 - r 2 (d6 2 + 1 



! ) (18) 



Here we have taken G = 
uses the symbols A 2 = e v 



= 1. Recall that one often 
g 00 and B 2 = e A = -g rr . 



The radial coordinate r is the luminosity distance and as 
before, the coordinate volume element is dV = Anr 2 dr. 
The proper volume element however is 



dV = B(r) dV = 4irr 2 B(r) dr 



(19) 



Note that since B(r) > 1, dV > dV, and this might be 
seen as the effect of "stretching" of space by gravity. The 
energy momentum tensor of the body in mixed tensor 
form is 



Ti = (p + p)u i u k +pgi 



(20) 



where u l is the fluid 4- velocity, p is the total mass energy 
density (excluding contribution due to self-gravitation), 
and p is the isotropic pressure. 

The total mass energy density of a cold fluid (excluding 
negative self-gravitational energy) is 



P = Pa + e 



(21) 



where pa = m^n is the proper rest mass energy density, 
m 7y is nucleon rest mass, n is nuclcon proper number 
density (not to be confused with polytropic index), e is 
the internal energy density, and g\ is the metric tensor. 
The explicit form of B(r) is [8] 

B(r) = [1 - 2M(r)/r}- 1 / 2 = (1 - 2m/r)- 1 ' 2 (22) 
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where 



M(r] 



r 

= m = 

Jo 



Awp r dr 



(23) 



The total energy of the system as perceived by a distant 
inertial observer, Soo, i.e., the gravitational mass of the 
fluid is 



l-R t-R 

M = I pdV = (p/B) dV 



o 



(24) 



o 



This total energy is also known as gravitational mass or 
Schwarzschild mass of the fluid although actually this 
should have been named by the name of Hilbert. Note 
that 

l-R [-R 

M = / (p/B) dV ^ pdV (25) 
Jo JO 

as one might have expected. The reason for this is that 
total mass energy includes not only local contributions 
from p but also the negative global contribution of self- 
gravitational energy. It is this negative latter contribu- 
tion which reduces the effective net proper energy density 
from p to p/B(r). 

On the other hand, the total proper energy content of 
the sphere, by excluding any negative self energy con- 
tribution , i.e., the energy obtained by merely adding 
individual local energy packets, is 



M proper = j p dV 




(26) 



The conventional definition GR self- gravitational en- 
ergy of the body is [7, 8] 

f R 

E G = M- M proper = (1 - V=g7r) P dV (27) 
Jo 

Since B(r) > 1 in the presence of mass energy, Eq is 
a -ve quantity, as is expected. Note however that the 
Eg defined by Eq.(27) is the sum of appropriate lo- 
cal (proper) quantities somewhat like the definition of 
M proper in Eq.(26); and is not defined with respect to the 
inertial observer Soo- Note also that, in contrast, gravi- 
tational mass M is indeed the mass-energy measured by 

Soo- 

The proper rest mass energy of the fluid is 



Po dV 



(28) 



If there are no initial anti-baryons or anti-leptons, Mo = 
tojv N, where N is the total number of baryons, and is 
therefore a conserved quantity. 

The proper internal energy of the fluid is 



U 



I 



e dV 



(29) 



As before, note that, neither Mq nor U is defined with 
respect to Soo- By using Eqs. (20-28), it can be verified 
that 



E G + U = M - M 



(30) 



From the above equation, it may appear that the GR 
equivalent of En, the binding energy, is 



E GR = M-M = E G + U 



(31) 



Suppose the gas was initially infinitely dispersed to in- 
finity. Further, let the gas molecules be rest with respect 
to Soo ■ Under such a case, the initial mass energy of the 
cloud is just the rest mass energy: 



M(t = 0) = M 



(32) 



And if the contraction of the cloud into a finite size would 
indeed release energy, one must have M — Mq < 0. But 
unlike the Newtonian noted by Tooper[9], it 

is not apparent from the definitions of Eq and U that 
Egr = Eg + U < 0! This is so because while in the 
Newtonian case, E^ and U are related through Eq.(7), 
there is no known relationship between Eq and U in the 
GR case. Further there might be examples when the 
occurrence of a negative Eqr "is not a sufficient condi- 
tion for instability of the system against an expansion 
to infinity" [9] . This means that Egr may not be the 
correct GR equivalent of a "binding energy" . Had Egr 
been the true GR equivalent of En , probably, one would 
have had equations similar to (8-9) involving Eq, U and 
Egr- But no such equations exist. Thus, although, intu- 
itively, one expects GR contraction to release radiation 
energy, one really cannot show it unlike the Newtonian 
case developed in previous section. 

In Eq.(30), we may note that, while, Eq, U and M 
are not defined with respect to (w.r.t.) Soo, M, on the 
other hand is defined w.r.t. Soo. And this may be the 
fundamental reason that Egr defined in Eq.(31) may not 
be the true "binding energy" of the fluid. 



ANOTHER DEFINITION OF FLUID MASS 

For any stationary gravitational field, total four mo- 
mentum of matter plus gravitational field is conserved 
and independent of the coordinate system used[10, 11]: 



/(T 



t0 + t l °) dV 



(33) 



where t lk is the energy momentum pseudo-tensor asso- 
ciated with the gravitational field. Further, the inertial 
mass (same as gravitational mass), i.e, the time compo- 
nent of the 4-momentum of any given body in GR can 
be expressed as[10-12] 



M 



/ m-Ti 

Jo 



-g d 3 x 



(34) 
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where 



g = -r 4 A(r)B(r) sin 2 9 



(35) 



is the determinant of the metric tensor go- and d 3 x = 
dr d9 d(f>. Since 

Tl = T 2 2 = T 3 3 = -p; T ° = p (36) 

it follows that [11, 12] 

M= (p + 3p)A(r)B(r)dV = / {p + ip)A(r)dV 
Jo Jo 

(37) 

When the body is bound and pb = Pb = for r > R, 
then, the foregoing integrals shrink to [11, 12] 

M= / (p + 'Ap)A(r)B(r)dV = / (p + 'Ap)A(r)dV 
Jo Jo 

(38) 

Although our result would not depend on splitting of 
the foregoing equation, (since we would simply equate 
the "total" expression of "inertial" and "gravitational" 
mass), we might nevertheless do so 

rR rR 

M = p A(r) dV+ ip A(r) dV (39) 
Jo Jo 

for the sake of obtaining physical insight. 



Global Definitions 

The Newtonian virial theorem (4) is essentially a 
statement of energy conservation involving negative self- 
gravitational energy and positive thermodynamic energy 
of the fluid. In the Newtonian case, there exists global 
inertial frames and a statement of energy conservation 
can be made in a trivial way. But in GR, even for this 
simplest case of a static fluid sphere, there is no global 
inertial frame. Thus the exercise of having global defini- 
tions of related energies and to enact their conservation 
is a highly non-trivial task. As is well known, for sta- 
tionary systems, however global energy can be defined 
in a meaningful way for asymptotically flat spacetimes. 
Further, when the energy is defined with reference to an 
observer at a spatial infinity (S^), we obtain the so-called 
ADM Mass[13]. Also global energy conservation can be 
meaningfully defined only with reference to Soo- 

Note that Eg and U are essentially summation over 
local appropriate values and not over the corresponding 
quantities measured by Soo- It may be mentioned that, 
if any locally measured energy is e, then the energy mea- 
sured by the far away inertial observer is the redshiftcd 
quantity 



Accordingly, the total mass energy content of the fluid, 
excluding any self-energy, as measured by an inertial 
frame such as the distant observer, 5^, is different from 



and is given by[8] (Eq. 11.1.19): 



M matter = jf * A{r) pdV = j A{r) B(r) p dV (41) 

Physically this means that local energy content in a given 
cell pdV is measured {redshifted) as A(r)pdV by the iner- 
tial observer Soo- 

In fact, the inertial observer Soo would see the rest 
mass energy, i.e., the proper energy, of a nuclcon too 
to be reduced by the same factor ^/goo- Hence, when 
the body is finite and not dispersed to infinity, the total 
rest mass energy of the body, as reckoned by the inertial 
observer 



M = / po A(r) dV 
Jo 



(42) 



Similarly, the redshifted global internal energy of the fluid 
as measured by is 



U 



f 

Jo 



e(r)A(r) dV 



(43) 



We can, now, quickly identify the 1st term on the RHS 
of Eq.(39) as M matter . Again bear in mind the fact that 
our eventual result would not depend on such identifica- 
tion or splitting because it would be obtained by equat- 
ing the "total" expressions for inertial and gravitational 
masses. 

Further, we see that, the 2nd term on the RHS of 
Eq.(39) is the global energy associated with pressure as 
perceived by Soo- Accordingly, we rewrite Eq.(39) as: 



M = M matter + M p 



(44) 



Mr) 



(40) 



Again recall that M too is defined only with respect to 
5qo. Having done this splitting, we are in a position to 
obtain the GR Virial Theorem, which is essentially an 
accounting of global energies involved in the problem. 
And since global energies, in GR, can be defined only 
w.r.t. the inertial observer Soo, ah relevant integrals must 
be defined w.r.t. the same observer. And this is what we 
have just done. 



GR HELMHOLTZ KELVIN PROCESS 

Let us simply transpose Eq.(39) (irrespective of its 
splitting) as 



/ p A(r) dV-M+ 3p A(r)dV = (45) 
Jo Jo 
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to reexpress as 



E n 



f 

Jo 



3p A(r) dV = 



Or, 



where 



E„ 



5oo 



dV = 



(46) 



(47) 



Eg — M — M mat ter 



I 



(AB - 1) p dV (48) 



Or else. 



E 



g - I (V-g009rr ~ 1) p dV 



(49) 



Since AB < 1 in the presence of mass-energy, we have 
Eg < as is expected. Clearly, we have, obtained, now 
an equation similar to (4) . It appears then that the above 
defined E g rather than the previously defined Eq is the 
true measure of self-gravitational energy as perceived by 
an inertial observer S^. This is so because Eq is not 
defined with respect to S^, the accountant for global 
energy. In contrast, both the components of E g , namely 
M and M matter arc defined w.r.t. Soo- 

Further, recall that, in GR, the effect of "gravitational 
potential" is conveyed by 500- But Eq (Eq.[27]) does 
not contain goo a t all. In contrast, E g indeed involves 
gravitational potential term goo = A 2 (Eq.[49]). 

What would be the value of true global GR self- 
gravitational energy E g in the Newtonian limit? 

To see this we consider a sphere with p = constant for 
which [8] 

A(r) = ^[3{1-2M/R)^ 2 - B(r)- 1 } (50) 

Using Eq.(22) in (50), we further see that 

A(r)B(r) = 1 [3(1 - 2M/R) 1 ' 2 B{r) - 1] (51) 

Again using Eq.(22) in (51), we obtain 

A(r)B(r) = i [3(1 - 2M/R) 1 ' 2 {1 - 2m/r)- 1 ' 2 - 1] (52) 

Now if we proceed to linearized gravity limit with 
M/R <C 1 and m/r <C 1, we will have 



A(r)B(r) = 1 - -(m/r - M/R)} (53) 

Using Eq.(52) in (48), we see that 
3 



E 9 = -o I (M/R-m/r)pdV 



(54) 



When we carry out this above integration with p = 
constant, we obtain 



E a 



-3M 2 

T~r 



(55) 



rpN 



Therefore, in the Newtonian limit, E g — Eq 
though in general E g and E g arc different. 

Now, using the thermodynamical relation (5) and 
Eq.(43) in Eq.(46), as before, we will have 



Eg + 3(7 ~1)U = 



(56) 



By direct comparison with Eqs.(4) and (7), we can easily 
identify Eqs.(48) and (56) as the appropriate GR version 
of static scalar virial theorem. Note that Eqs.(46) and 
(56) naturally reduce to their Newtonian forms, Eqs.(l) 
and (7) for sufficiently weak gravity with goo ~ —Qrr ~ 1- 
Thus we may interpret the existence of the Newtonian 
virial theorem too as due to equivalence of "gravitational 
mass" and "inertial mass" . 
From Eq.(56), we obtain 



U = 



1 E 
3(7-1) 9 



(57) 



If the fluid undergoes quasistatic contraction and A de- 
notes the associated changes in relevant quantities, then 
we will have 



Af> = —^—AE Q = , +l - 

3(7-1) 9 3(7-1) 



AE n 



(58) 



Here we have used the fact that since E g < it must 
decrease for contraction. Thus as is expected, the in- 
ternal energy of the fluid must increase for gravita- 
tional contraction. If the appropriately averaged value 
of A(r) = goo during this contraction, the increase in 
proper internal energy would be 



AU 



AU 

900 



\AE n 



3(7 - 1) 5oo 



> 



(59) 



Since goo < 1 in the presence of gravity, this means that, 
the rate of increase in proper internal energy would be 
higher than in the corresponding Newtonian case. 

As we look back at Eq.(58), the increase in the value 
of I .E g |, namely |AE g | is not fully accounted for by the 
increase in the value of U. Note that in the absence of 
initial antiparticlcs, the contribution of rest mass-energy 
is unaffected during the process. Therefore, for the sake 
of global energy conservation, as reckoned by the inertial 
observer S^, the fluid must radiate out an amount of 
energy +AQ given by 



AQ 



1 



1 



3(7-1) 



\AE n 



(60) 
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in order to be able to contract. Note, as before, that 
7 > 4/3 in a strict sense, as long as particle momenta are 
finite. To see that in the GR context too, that 7 = 4/3 
implies singular situation, look at the 1st row of Table I 
of [14] which shows that in this case again R max = 0, 
where R ma x is the maximum possible radius of the con- 
figuration. Further, for 7 = 4/3, the next entry in the 
same row shows that R = where Rq = 2GM/c 2 is 
the Schwarzschild radius. This implies that for 7 = 4/3, 
one has Rq = 0. This latter fact implies that even in 
GR, the total mass energy Mc 2 = just as En = for 
7 = 4/3 (Eq.[9]). Occurrence of R ma x = implies a fluid 
sphere that has collapsed to a singular point. And as 
per Ref.(14), the singular configuration then would have 
zero mass energy. Incidentally, in the classic paper [13], 
Arnowitt, Dcser and Misncr too found that a neutral 
"point particle" has zero "clothed mass" . Then Chan- 
drasekhar's exercise[14], in addition, suggests that if the 
fluid would attain such a singular state, the value of 
7 — > 4/3. This is also in perfect agreement with the 
notion that a singular state should be infinitely hot with 
complete domination of radiation energy over rest mass 
energy [6]. 

In fact, in one would misconceive of a situation with 
7 < 4/3, irrespective of whether it is a Newtonian or a 
GR case, one would have to ensure injection of energy 
into the system to let it collapse. This would mean that 
in the absence of external injection of energy, the system 
would not contract /evolve at all in defiance of basic tenet 
of global gravitation. 

Thus, one must indeed have 7 > 4/3 and, in a strict 
sense, there cannot be any adiabatic gravitational contrac- 
tion. Consequently, all strictly adiabatic gravitational 
collapse studies are of only academic interest. Thermo- 
dynamically, the global specific heat of the contracting 
fluid is negative because dQ < while the temperature 
and internal energy of the fluid increase. 

GENERAL RELATIVISTIC BINDING ENERGY 

Eq.(60) suggests that, we might isolate a quantity 

as the total energy of the system excluding any rest mass 
contribution. This is actually the "Binding Energy" of 
the gravitating system, defined by a given distribution of 
mass-energy. 

By using Eq.(57), it is seen that 

E = E g + U (62) 

Further using Eqs.(39), (42), (43) and (48), it also tran- 
spires that 

E = E CJ + U = M — Mq (63) 



We may see that, unlike in Eq.(30), all the quantities 
involved in Eq.(63) are defined w.r.t. and which sug- 
gests that E is indeed the true binding energy of the 
fluid. The reader is again requested here to appreciate 
the subtle point why the true binding energy of the fluid 
is given by Eq.(63) rather than by Eq.(31).: 

Though Mq is a conserved quantity, once the fluid is 
contracted into a finite size, it gets dissociated from the 
inertial frame 5^. And S^, who is doing the global en- 
ergy accounting, sees the locally defined rest mass energy 
to be rcdshiftcd to Mq rather than as Mq. On the other 
hand, the total mass-energy of the fluid, again defined 
w.r.t. Soo is M. Therefore, the global binding energy 
of the fluid, i.e., the difference between the total mass 
energy and rest mass energy as seen by the same inertial 
observer is M — Mq. Note that the existence of equa- 
tions (62-63) does not depend on such interpretations 
because they, in any case, crept up spontaneously. 

Although, in a Newtonian case, the notion of a bind- 
ing energy always existed, to the best of our knowledge, 
such a notion was never before properly derived in the 
GR context. A relativistic "bound system" may thus be 
defined as one having E < 0, and an "unbound system" 
will have E > 0. 

While, in the Newtonian case, "Total Energy" is the 
binding energy En, in GR, total global energy, as mea- 
sured by Soo, always is E = Mc 2 . 

DISCUSSIONS 

The important idea of Hclmholtz and Kelvin, devel- 
oped in the 19th century, that gravitational contraction 
should both raise the internal energy and cause the fluid 
to radiate was always expected to be valid irrespective of 
the strength of the gravity. However, the original deriva- 
tion to this effect was made in the framework of extremely 
weak, i.e., Newtonian gravity. We showed here that even 
for arbitrary, strong gravity, this process indeed remains 
valid. In fact, as shown by Eq.(59), the process becomes 
even more effective compared to the Newtonian case as 
gravity increases. Pictorially, we may think that stronger 
gravity churns out more radiation from matter even in 
the absence of chemical or thermonuclear energy sources. 
A similar conclusion is supported by a recent work which 
shows that the ratio of radiation energy density to rest 
mass energy density of a self-luminous contracting ob- 
ject is proportional to its surface redshift z [6]. When 
z<l, the self-gravitating contracting object is "matter 
dominated", i.e., po 3> p r , but when, z 3> 1, the object 
becomes radiation dominated : p r ^> Po like the very 
early Universe [6]. The present study provided an addi- 
tional explanation for this result. These studies show 
that the actual fate of radiative physical gravitational 
collapse could be radically different from traditional pic- 
tures of continued gravitational collapse inspired by the 
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pressureless dust collapse where a Black Hole (BH) or a 
Naked Singularity is catastrophically formed in a finite 
comoving proper time. Traditional GR collapse studies 
are usually done by (i) assuming dust models with p = 
even when the fluid is supposed to have collapsed to sin- 
gularity or (ii) considering pressure but neglecting all 
heat transport, dQ = 0. But as shown by Eqs. (58) and 
(60), dQ = implies (a) dU = and (b) dE g = 0. The 
condition (a) is satisfied only for dust and thus despite a 
formal consideration of existence of pressure, in the con- 
text of collapse, a fluid satisfying condition (a) becomes 
similar to a pressureless, internal energy less dust. The 
condition (b) is not satisfied even for a dust unless it 
has M = U = Eg = fixed = at the beginning of the 
collapse. But no isolated fluid with finite size can have 
M = (the Universe may, however, have M = even 
being of finite or infinite extent). In several numerical 
studies of supposed radiative collapse, one implicitly or 
explicitly assumes Q -C Mqc 2 . At the advanced stage 
of collapse, this assumption fails [6] and such cases effec- 
tively become similar to case (ii) of adiabatic continued 
collapse valid for M = 0. 

In contrast, in a breakthrough research on physical 
gravitational collapse, Herrera and Santos[l] have shown 
that the force exerted on the collapsing fluid by the out- 
ward propagating heat / radiation may stall the continued 
collapse process and formation of either a finite mass 
BH or Naked Singularity may be averted. Herrera, Di 
Prisco and Barreto[l] have successfully made a numeri- 
cal model of continued collapse to substantiate this path- 
breaking idea. Such ideas are consistent with the model 
independent generic studies [6] which show that contin- 
ued catastrophic collapse indeed degenerates into a ra- 
diation pressure supported hot quasistatic state called 
eternally collapsing objects because of outward force due 
to collapse generated radiation at extremely deep grav- 
itational potential wells, z ^> 1 where z is the surface 
gravitational redshift of the collapsing object. It is be- 
cause of the resultant reduction in the value of M due to 
continuous radiation outpour that no apparent horizon 
or event horizon is formed until M = R = 0[6]. 

Finally, the GR definition of "binding energy" of a 
static fluid is M - M rather than M - M Q . 

Newtonian HK process is a direct sequel of static New- 
tonian virial theorem. Similarly, we needed to derive 
the exact GR version of the static virial theorem. This 
GR virial theorem involved globally defined quantities 
measured w.r.t. the same inertial observer S^. In gen- 
eral, the notion of global energy is far from transparent 
and unique in GR. For example, for non-static and non- 
sphcrically symmetric systems or charged systems, there 
could be various notions of "energy" and "mass" ; to ap- 
preciate this one may have a look at a recent long review 
paper [15]. However the present paper must not be con- 
fused with such studies. The aim of this paper was not 
at all to define any new definition of either "mass" or 



"global energy" from any preferred theoretical perspec- 
tive, correct or incorrect. This is so because, unlike a 
generic case, the definition of "global mass energy" of 
a chargeless static spherically symmetric fluid (measured 
by Soo) is very well known since long[10, 12]. And we just 
appealed to the Principle of Equivalence to equate the 
already well known expressions for "gravitational mass" 
(Schwarzshild mass) and "inertial mass" - the time com- 
ponent of linear 4-momentum. It is this simple opera- 
tion which yielded the GR virial theorem and GR HK 
mechanism (Eqns. [12-17]). It is the same principle of 
equivalence which demanded that, from energy conser- 
vation considerations, all the relevant globally defined 
energies are defined w.r.t. the unique inertial frame 
rather than w.r.t. a series of (infinite) proper frames. To 
the best of our knowledge, this is the maiden derivation 
of HK process using GR. And this is also the maiden 
proper GR explanation for the intuitive notion that a 
sclf-gravitating object has effective global "negative spe- 
cific heat" in Einstein gravity too. 

For further appreciation of the "global quantities" in- 
volving y/goo in our study which arose spontaneously and 
not imposed from new theoretical perspective, we recall 
that the "Poisson's Equation" in GR has the form[16] 

V 2 V^ = 47rGV5^(p + 3p) (64) 

And only when one moves to weak gravity with 
GM/r <C 1 and ^/<?oo w 1 — GM/r w 1, one obtains 
the more familiar form 

V 2 = 4ttG (p + 3p) (65) 

where the weak "gravitational potential <f> ~ GM/r. Of 
course, Eq.(64) would also degenerate into Eq.(65) in a 
local free falling frame where goo = 1 ■ But for fluids hav- 
ing finite pressure there cannot be any such global free 
falling frame and therefore Eq.(65), in the present con- 
text, can be recovered only for weak gravity (in any case 
we are dealing with a static fluid. Formally a strict adi- 
abatic collapse is possible only for a pressureless "dust" 
with p = U = 0, though, physically, M = in such a 
case. 

Since virial theorem is important for the study of com- 
pact objects and gravitational contraction, the exact rel- 
ativists virial theorem obtained here could be useful for 
relativistic astrophysics, either now or in future. 
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